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sensitive to small variations in interatomic distances within the coordination shell of the absorbing 
atom. The experiments were carried out on single crystals of ruby and sapphires. Linear dichroic 
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basis set, and the continued fraction for the absorption cross section. 
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I. INTRODUCTION 

The presence of impurities in a crystal can influence its 
mechanical, magnetic, semiconducting, superconducting, 
dielectric, spectroscopic and transport properties. To cite 
just a few specific examples, impurities can improve the 
voltage holdoff capability of insulating materials^, they 
are critical for the optical properties of most laser mate- 
rials and gemstones, and they turned out to be the secret 
of ancient Damascus steel bladesjS. 

Still, the physics of impurity systems is not well de- 
veloped. The local crystallographic structure around im- 
purities is unknown in most cases, although it is an es- 
sential piece of information to undestand their influence 
on the physical properties of the host and to carry out 
ab initio simulations of these materials. This situation is 
due to the experimental and theoretical difficulties that 
are met to obtain microscopic information on impurity 
systems. From the theoretical point of view, quantum 
calculations of impurity systems require the use of large 
supercells that have long not been computationally man- 
ageable. From the experimental point of view, impurity 
systems can only be measured with methods which can 
selectively probe certain atomic species. Among such 
methods, many of them, such as electronic paramagnetic 
resonance, give only indirect information on the position 
and nature of the atoms surrounding the impurity. Due 
to the advent of third-generation synchrotron radiation 
facilities, x-ray absorption spectroscopy (XAS) is now 
able to investigate impurities in solids. 

The extended x-ray absorption fine structure (EX- 
AFS) region of a XAS spectrum provides quantitative 



information about the short range organization around 
the absorbing atom (coordination number, interatomic 
distances). The x-ray absorption near-edge structure 
(XANES) region usually gives qualitative information 
about the atomic arrangement up to the medium range 
order. XANES is sensitive to the electronic structure 
(bonding) and probes the empty states of solids. In the 
case of dilute elements in single crystals, especially at 
the K edge of 3d transition elements, collecting EXAFS 
spectra with a good signal-to-noise ratio is a difficult task, 
because the intensity of diffraction peaks becomes much 
larger than the EXAFS signal. Taking the example of 
impurities in corundum a-Al203,^ good EXAFS spectra 
could be collected for 10000 wt.ppm Cr'^+ impurities, but 
the usable energy range obtained for 1500 wt.ppm Fe'^+ 
impurities was already quite limited, and the concentra- 
tion of coloring impurities is often an order of magnitude 
lower than this. In such cases, XANES becomes a prac- 
ticle and precious technique to get information about the 
local structure of the absorbing atom,i However, the an- 
alyzis of the experimental XANES data is not straight- 
forward because of the photoelectron multiple-scattering 
processes that occur in the near-edge region. Ab initio 
XANES simulations are then required to relate the ex- 
perimental spectral features to the local geometry around 
the aborbing atom. 

Recently, an efficient first-principle approach based on 
plane-wave pseudopotential formalism has been devel- 
oped to calculate K-edge XANES spectra.^ Here we ap- 
ply this method to the K-edge of substitutional param- 
agnetic impurities in aluminum oxide. More precisely, we 
investigate the structural modifications of the corundum 
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crystal structure induced by the presence of substitu- 
tional Ti'^+, Cr^+ or Fe'^"'". This system is chosen because 
corundum containing transition metal impurities is im- 
portant for laser applications and as gemstones. Colorless 
corundum (a-Al203) becomes red ruby, pink sapphire, or 
yellow sapphire when a small amount of Cr'^+j Ti'^^, or 
Fe'^^ ions substitute for Al'^^ ions, respectively.^-^ The 
color of blue sapphire is due to the presence of (Fe-Ti) 
pairs in a-Al203. 

For this investigation, polarized XANES spectra are 
measured and calculated. Corundum belongs to the trig- 
onal Bravais lattice, the optical axis being parallel to the 
[111] direction of the trigonal unit cell (i.e. to the c-axis 
of the hexagonal cell) . Corundum is then a dichroic com- 
pound in the electric dipole approximation. The dichroic 
signal is a direct signature of the departure of the im- 
purity site from octahedral symmetry.-^ XANES calcu- 
lations are performed from structural models resulting 
from ab initio energy minimization calculations. The 
agreement of our experimental and theoretical spectra 
demonstrates that precise structural information can be 
extracted from the anglular-dependent XANES spectra, 
provided the electronic potential is accurately modelled. 



II. EXPERIMENTS 

A. Crystallographic structure of Q-AI2O3 and 
samples description 

Corundum (q!-A1203) belongs to the R3c (D^^) space 
groupi^ The trigonal unit cell contains two AI2O3 formula 
units. The atomic environment of aluminum is illustrated 
in Fig. n The atomic site of aluminum is a distorded oc- 
tahedron with a 3 (C3) local point symmetry. The AlOg 
octahedron is characterized by two different Al-0 inter- 
atomic distances. If we consider that the absorbing atom 
is denoted by Al in Fig.^ the nearest three oxygen atoms 
are labeled Oi in Fig.Q] The farther three oxygen atoms, 
labeled O2, form a face shared by two octahedra along 
the three-fold symmetry axis. Beyond the coordination 
shell of aluminum, the next two neighbors are aluminum 
atoms, labeled Ali and AI2 in Fig. and are relative 
to face-shared octahedra and edge-shared octahedra, re- 
spectively. 

Three synthetic single crystals of doped a-Al203 (ob- 
tained by the Verneuil process) are used for this study: a 
red ruby (a-Al203:Cr^+), a blue sapphire (a-Al203:Fe- 
Ti) and a pink sapphire (Q!-Al203:Ti'^+). The impurity 
concentration in each sample is given in tableQ] They are 
measured with the Cameca Microbeam electron micro- 
probe at the CAMPARIS analytical facility of Paris. A 
30 kV acceleration with a 15 uA beam current, defocused 
to lOiim, is used. X-ray intensities are corrected for dead- 
time, background, and matrix effects using the Cameca 
ZAF routine. The standards used are q;-A1203, a-Cr203, 
Q;-Fe2 03 and Ti02. The blue sapphire sample is cho- 
sen to represent the model compound of iron impurity 
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FIG. 1: Aluminum atomic site in corundum. Left: the corun- 
dum structure is represented as a stacking of AlOe octahedra 
(the trigonal cell and the three-fold symmetry axis are indi- 
cated). Right: zoom of the black- framed region of the left 
panel. The aluminum site with its Oi, O2, Ali and AI2 first 
four neighbors are indicated. 

in corundum (Q!-Al203:Fe'^~''). It is indeed preferred to a 
yellow sapphire sample, that only contains iron impuri- 
ties, because synthetic yellow sapphires are known to be 
inhomogeneousiS Consequently, in this study, we assume 
that the low concentration of titanium (750 wt.ppm) 
compared to that of iron (1500 wt.ppm) does not affect 
the site relaxation around iron. 



TABLE I: Impurity concentration (in wt.ppm) in the three 
doped a-Al2 03 samples. 



Samples 


Cr (wt.ppm) 


Fe (wt.ppm) 


Ti (wt.ppm) 


Red ruby 


10000 (± 1000) 






Blue sapphire 




1500 (± 50) 


750 (± 30) 


Pink sapphire 






540 (± 30) 



The samples are cylindrically shaped: 15 mm diame- 
ter and 3 mm thickness for the ruby sample, and 5 mm 
diameter and 1 mm thickness for the sapphire samples. 
They are cut so that the optical axis is in the plane of 
the disk surface, and the [101] direction of the trigonal 
cell is orthogonal to the disk surface. 



B. X-ray absorption measurements 

X-ray absorption measurements are carried out at the 
European Synchrotron Radiation Facility (ESRF) on the 
ID12 beam line, which is dedicated to polarization depen- 
dent spectroscopies. ^^'-J-' The monochromatic x-ray beam 
is obtained through a double Si(lll) crystal monochro- 
mator, with a resolution AE/E ~ 10^**. The oriented 
samples are placed on a rotating holder with the normal 
of the disk surface parallel to the x-ray wave vector and 
to the rotation axis. The fluorescence intensity is mea- 
sured by eight detectors symmetrically set around the 
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x-ray beam. 

XANES spectra are recorded for several angles 9 be- 
tween the optical axis of the sample and the x-ray polar- 
ization vector. We note (T|[ and the absorption cross 
section obtained with = 0° and = 90°, respectively. 
In the electric dipole approximation, the isotropic sig- 
nal is given by (cr|| -|- 2a_L)/3, and the dichroic signal is 
given by the difference a\\ — a^. The measurement of 
x-ray absorption spectra of impurities in crystals is im- 
peded by the presence of large diffraction peaks. In order 
to decrease the intensity of diffraction peaks and elastic 
scattering, V and Mn filters (50 /zm width) are used for 
the Cr K edge in ruby and the Fe K edge in blue sap- 
phire, respectively. By so doing, the spectra recorded for 
the two angles = 0° and — 90° are not affected by 
diffraction peaks. On the contrary, no filter is available 
for the Ti K edge and the absorption was measured for 
100 angles from to 360°. The a\\ and spectra are 
reconstructed from this amount of data by using a Fourier 
decomposition to remove the contribution of diffraction 
peaks to the spectra. 



C. Calculation methods 

The starting point of any XANES ah initio calculation 
is a structural model. Two structural models are consid- 
ered for each compound. The first one, called non-relaxed 
model, corresponds to the substitution of one aluminum 
atom by one impurity Ti, Cr or Fe in the q;-A12 03, with- 
out relaxing its structure. The second one, called relaxed 
model, is the result of ah initio energy minimization cal- 
culation carried out on the non-relaxed model. This cal- 
culation is performed using Car-Parrinello molecular dy- 
namics with simulated annealing^^ as implemented in the 
CPMD code.^^ The XANES calculations are then carried 
out within the scheme of the PARATEC coder^ using the 
formalism described in Refj^ Both codes are based on the 
density functional theory (DFT) within the spin polar- 
ized local density approximation (LSDA) , and use plane- 
wave basis set and norm-conserving TrouUier-Martins 
pseudopotentialsi^ in the Kleinman-Bylander form>^ In 
the following, we give all the parameters used in both 
kinds of calculation. 



1. Structural model construction 

The structural models are constructed from the lattice 
parameters and atomic positions resulting from a previ- 
ous ah initio calculation, performed by Duan et al}'^'^^ in 
a-Al2 03. In Refiii, the trigonal unit cell parameters are 
an = 5.11 A and a = 55.41°, the aluminum atoms are 
in 4 c symmetry sitesiS with x = 0.352 and the oxygen 
atoms are in 6 e symmetry sitesiS. with x = 0.555. We 
verify by using CPMD that this corundum optimized struc- 
ture obtained by Duan et al. corresponds to atomic forces 
less than 2.10-^ Ry/A. Since we use periodic bound- 



ary conditions, we consider 2x2x2 supercells. The 
supercells contain 80 atoms: 48 oxygen atoms, 31 alu- 
minum atoms and 1 transition metal atom (Ti, Cr or Fe) 
in substitution for aluminum. With such large super- 
cells, the interaction between two impurities (belonging 
to two neighboring cells) is negligible. These supercells 
define the non-relaxed models mentioned above. The re- 
laxed models are obtained by minimization of the energy 
of the non-relaxed supercells. All atoms are allowed to 
relax, while the lattice constants are fixed. The impuri- 
ties considered in this study are known to be in a high- 
spin state, ^4 therefore the spin multiplet imposed on the 
trivalent impurities are S=i for T\^+ ([Ar]4s°3di), S=| 
for Cr3+ ([Ar]4s°3d3), and S=f for {[kr]As°2,d^). 
The parametrization of the norm-conserving TrouUiers- 
Martins pseudopotentials is given in Tabled The wave 
functions and the charge density were expanded in plane 
waves with cutoff of 80 Ry and 320 Ry, respectively. Since 
the supercell is rather large, and since the systems are in- 
sulating materials, the Brillouin zone is only sampled at 
the r point. 

TABLE II: Parametrization used for the generation of the 
pseudopotentials. The core radii of the valence states are 
indicated between parenthesis in A. 

atom Al O Ti Cr Fe 

valence 3s (1.06) 2s (0.77) 3s (0.58) 3s (0.53) 3s (0.48) 
states 3p (1.06) 2p (0.77) 3p (0.90) 3p (0.79) 3p (0.90) 

3rf (1.06) 3d (0.90) 3d (0.79) 3d (0.90) 

local part d p d d d 



2. XANES calculations 

The method used for XANES calculations has been al- 
ready described in Rcf.^-^^ Therefore its main aspects 
are recalled here. The method uses periodic bound- 
ary conditions, plane-wave basis-set, pseudopotentials 
and reconstructs all-electron wave functions within the 
projector augmented wave (PAW) framework. In or- 
der to allow the treatment of large supercells (hundreds 
of atoms), the scheme uses a recursion method to con- 
struct a Lanczos basis and then compute the cross sec- 
tion as a continued fraction^SiSi in the electric dipole 
and quadrupole approximationsiS^ Electric quadrupole 
transitions are relevant for the pre-edge region, in par- 
ticular at the i^-edge of transition metals. The absorp- 
tion cross section is calculated beyond the "muffin-tin" 
approximation, that is known to limit the applications 
of the multiple scattering theory traditionally used for 
XANES simulations.2SiaL2a 

XANES spectra are computed at the transition metal 
K edge for both relaxed and non-relaxed models of the 
three compounds, using the 80- atom supercells described 
in the previous subsection. For the pseudopotentials con- 
struction, we again use the parametrization given in Ta- 
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FIG. 2: Energy dependence of the 7 parameters used for the 
convolution of the XANES spectra at the Ti, Cr and Fe K 
edges in doped a-Al2 03. The zero in the energy scale refers 
to the highest occupied level in the corresponding structure. 



blem In order to take into account the core- hole effects 
in the calculations, the Ti, Cr and Fe pseudopotentials 
are generated with only one Is electron. The spin multi- 
plet degeneracy is set by imposing the number of up and 
down states given by the CPMD code. Convergence of the 
XANES theoretical spectra is reached for the following 
set of parameters: a 64 Ry energy cutoff and one fc-point 
for the self-consistent spin-polarized charge density calcu- 
lation, and 32 fc-points for the absorption cross section 
calculation. An energy dependent broadening parameter 
7 is used in the continued fraction (see Refi^^Si) in or- 
der to account for the two main photoelectron damping 
modes. '^'^ This width 7 corresponds through the time- 
energy uncertainty relations to a lifetime of the photo- 
electron. At low energy, the lifetime of the photoelectron 
is only limited by the core-hole lifetime within the plas- 
mon pole approximation. The parameter 7 was then set 
to a constant, which corresponds to K level full width at 
half maximum given by Krause and Oliver At higher 
energy, the kinetic energy of the photoelectron is great 
enough so that it can interact with the electron gas of 
the system. Its amplitude is damped, provoking a sharp 
decrease in the lifetime. Therefore the 7 parameter was 
set to rapidly increase, and then to have a smooth linear 
behaviour. The 7 energy-dependent parameters used in 
this study are displayed in Fig. |2| 



III. 



RESULTS AND DISCUSSION 



This section is organized as follows. First, we present 
the results of the structural relaxation around the substi- 
tutional Ti3+, Cr3+ and Ye^+ ions in a- AI2O3 (SeclmSl. 
Second, we discuss the influence of the core-hole ef- 
fects in the case of the Cr K-edge in ruby fSec lIIIBjl . 
Third, we confront with experiments the calculated spec- 
tra obtained with both relaxed and non-relaxed models 

rsec inm . 



A. Structural model analysis 

The relaxed models of a-Al2 03:Cr^+ and a- 
Al203:Fe^+ are already described and discussed in a pre- 
vious work."^ However, in order to compare with the re- 
laxed model of a-Al203:Ti^+, most of the results of Refi^ 
are recalled here. 

The analysis of the relaxation can only be done by care- 
ful comparison of the relaxed and non-relaxed supercells. 
This is achieved by comparing clusters that are built 
from the supercells using the method detailed in Refj^ 
These clusters are centered on the impurity and contain 
66 atoms (5.2 A radius). In order to avoid the influence 
of the impurity displacement in the relaxed model, we 
define the mass center D, of each cluster by calculating 
OVt = ^ Si ^iOMi, where O is the impurity posi- 
tion, Ali refers to all the atoms within a given cluster 
except the impurity and mi is the mass of the atom i. 
For each atom i around the impurity, we evaluate the 
norm of the displacement vector, Vi, 
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where 

nates of fiMi vector in the non-relaxed and relaxed 
clusters, respectively. The angular relaxation was 
also determined. It was defined by the quantity 

— ^non-relaxed ^relaxed: where O^elaxed and rion-relaxed 

are the angles between the fiM^ directions and the C3 
axis in the relaxed and non-relaxed clusters, respectively. 
The comparison of the relaxed and non-relaxed clusters 
leads to the following conclusions. 
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FIG. 3: Angular relaxation around the paramagnetic impuri- 
ties in doped a-Al2 03: 59 angle (see text for explanation) as 
a function of distance from the impurity in the relaxed model. 

First, the impurity site symmetry is conserved: the Ti, 
Cr and Fe atoms are still in a C3 point symmetry site, 
characterized by two kinds of M-Oi and M-O2 distances 
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TABLE III: M-0 and M-Al bond lengths in Q:-Al2 03:M^+ (M=Ti, Cr, Fe) issued from the energy minimization calculations 
presented in Sec lIII A1 These bond lengths are compared with the experimental and calculated Al-0 and Al-Al interatomic 
distances in Qf-Al2 03 and to the M-O and M-M distances in q:-M2 03. It is worth noting that the distances indicated for 
C1-AI2O3 relative to Ref;lA are the M-O and M-Al distances in the non-relaxed models. 

Compound Interatomic distances (in A) 

cation-oxygen cation-cation 

a-Al203 (expiS) 
a-Al2 03 (calcii^) 
a-Al203:Ti^+ (this calc.) 
a-Al203:Cr3+ (calci^) 
a-Al203:Fe^+ (calci^) 
a-Ti203 (expiS) 
a-Cr2 03 (exp4a) 
Q-Fe203 (exp4&) 



dAl-Oi 


=1.86 


dAi-02 


=1.97 


dAi- All =2.65 


dAi-Ai2=2.79 


dAl-Oi 


=1.84 


dAi-02 


=1.95 


dAi-Aii=2.62 


dAi-Ai2=2.77 


dxi-Oi 


=1.97 


dTi-Oa 


=2.07 


dTi-Aii=2.72 


dTi-Ai2=2.83 


dcr-Oi 


=1.95 


dcr-Oa 


=2.00 


dcr-Ali=2.64 


dcr-Al2=2.80 


dpc-Oi 


=1.92 


dFo-02 


=2.10 


dFc-Aii=2.72 


dFG-Ai2=2.83 


dxi-Oi 


=2.01 


dTi-02 


=2.08 


dTi-Tii=2.59 


dTi-Ti2=2.99 


dcr-Oi 


=1.96 


dcr-02 


=2.01 


dcr-Cri=2.65 


dcr-Cr2=2.88 


dpc-Oi 


=1.94 


dFe-02 


=2.11 


dFe-Fei=2.90 


dFe-Fc2=2.97 



(M — Ti, Cr or Fe). The angular relaxation is found 
to be very small (i.e. less than 0.8 degrees). This is 
illustrated in Fig. O which gives 56 as a function of the 
distance between the impurity and its neigboring atoms 
in the relaxed clusters. 

Second, a displacement of the impurity from the initial 
aluminum site is observed in the three relaxed models. 
This displacement occurs along the C3 axis, but its ab- 
solute value depends on the impurity. Chromium is dis- 
placed by 0.03 A towards Ali atom, titanium is displaced 
by again 0.03 A but in the opposite direction, and iron 
is further moved (0.09 A towards Ali atom). The small 
displacement of chromium is in agreement with various 
experimental studies of the literature i^24Si2^i2Si2& The re- 
sult for iron is consistent with an analysis of electronic 
paramagnetic resonance (EPR) experiments, that gives 
a displacement of 0.04 ± 0.02 A,^ 
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FIG. 4: Radial relaxation around the paramagnetic impuri- 
ties in doped a-Al203: norm of atomic displacement vectors 
(oxygen atoms in solid line with circles and aluminum atoms 
in dashed line with triangles) as a function of distance from 
the impurity in the relaxed model. 

Third, an increase of distances within the coordina- 
tion shell of the impurity is observed. This result is 



expected since the ionic radii of Ti, Cr and Fe in oc- 
tahedral siteSa (rTi3+ = 0.670 A, rc,.3+ = 0.615 A and 
rp(,3+ = 0.645 A) are all greater than the ionic radius 
of aluminum (rAi3+ = 0.535 A). The calculated M-Oi 
and M-O2 bond lengths in doped Q-AI2O3 were found to 
be close to M-Oi and M-O2 bond lengths in the corre- 
sponding metal oxides, q;-M203 (Q;-Ti203, a-Cr203 and 
a-Fe203). On the contrary, the distances M-Ali and M- 
AI2 were less affected by the relaxation process: they are 
closer to Al-Ali and AI-AI2 bond lengths in corundum 
than the M-Mi and M-M2 distances in the corresponding 
oxides a-M203. Tablmlgathers the M-Oi, M-O2, M-Ali 
and M-AI2 distances in q;-A12 03:M^+ deduced from the 
calculation, and compares them with the corresponding 
distances in Q!-Al203, a-Ti203, Q!-Cr203 and Q;-Fe203, 
resulting from x-ray diffraction experiments. One can re- 
mark in Tab.Hnlthat the Ti-Oi and Ti-02 bond lengths 
in Q;-Al203:Ti''^+ (1.97 A and 2.07 A) are quite similar to 
those in Q;-Ti203 (2.01 A and 2.07 A). Besides, the Ti- 
Ali and Ti-Al2 distances in a-Al203:Ti3+ (2.72 A and 
2.83 A) are closer to the Al-Ali and AI-AI2 distances 
in corundum (2.65 A and 2.79 A) than to the Ti-Tii 
and Ti-Ti2 distances in a-Ti203 (2.59 A and 2.99 A). 
One can notice that the experimental distances Cr-Cri 
in Q;-Cr203 and Al-Ali in q;-A1203 are both equal to 
2.64 A. Thus one expects that Cr-Ali distance in ruby 
(a-Al203:Cr3+) should be around 2.64 A. The calculated 
distances of Q;-Al203:Cr^+are slightly larger (-^0.015 A) 
than those determined from EXAFS measurements on 
powder samples^^ and in good agreement (-1-0.005 A) 
with EXAFS measurements performed and on the same 
single crystal as the one used in the present study^ 
They are also consistent with data extracted from opti- 
cal spectroscopy!^ For a-Al203:Fe^"'", the calculated dis- 
tances are in agreement with EXAFS experiments^ and 
with EPR resultsi^ 

Finally, the M-O and M-Al distances deduced from 
the ab initio relaxation calculations suggest that the 
structural relaxation around the transition metal impu- 
rities is quite local and mainly concern the coordina- 
tion shell. This point is confirmed by the calculations of 



6 



the atomic displacements between the relaxed and non- 
relaxed. Fig. 21 represents the norm \Vi\ as a function of 
the distance between the impurity and its neighbors in 
the relaxed model. It clearly shows that the oxygen co- 
ordination shell absorbs almost the relaxation. The dis- 
placement of atoms farther than 2.5A from the impurity 
is less than 0.035 A. Comparing the three compounds, 
one also observes that the relaxation is slightly larger 
around titanium and iron than around chromium. 



B. Core- hole effects 

The influence of the core-hole can be described by cal- 
culating a self-consistent potential around a core-hole or 
by solving a Bethe-Salpeter equation. ^"•'^^ For impurity 
systems, solving a Bethe-Salpeter equation is probably 
beyond the reach of available computer programs. More- 
over, it was already observed that both methods give 
quite similar resultsi^ Therefore, we shall describe the 
influence of the core-hole by calculating a self-consistent 
potential in the presence of a Is core-hole. 

Calculated spectra performed with or without the pres- 
ence of a Is core- hole on the chromium absorbing atom 
are compared in Fig. It clearly appears that the core- 
hole is needed to obtain a good agreement with the ex- 
perimental curve. In particular, the intensity and the 
energy positions of the features within the 6005-6017 eV 
energy range are not well reproduced if the core-hole is 
not taken into account in the calculation. We have no- 
ticed similar behaviors due to the presence of the core- 
hole at the titanium and iron K-edges in Q;-Al203:Ti''+ 
and a-Al2 03:Fe'^+, respectively. 

On the other hand, the presence of the core-hole has a 
non-negligible impact in the pre-edge region. This is illus- 
trated in the inset of Fig.|5l which represents the prc-cdgc 
region (5985-5998 eV) of the Cr K edge in a-Al203:Cr3+. 
An interpretation of this region in terms of group theory 
has been detailed in Ref.''^ In the pre-edge region both 
electric dipole (Is ^ p) and quadrupole (Is — > 3d) tran- 
sitions occur. Since the site of the impurity in doped a- 
AI2O3 is not centrosymmetric, the p states are present in 
the pre-edge through p ^ d hybridization, electric dipole 
(Is — > p) and quadrupole (Is 3d) transitions occur. 
Since the site of the impurity in doped a-Al2 03 is not 
centrosymmetric, the p states are present in the pre-edge 
through p—d hybridization. The calculated curves shown 
in the inset of Fig. are the sum of both contributions. 
For the calculation of electric quadrupole transitions, the 
direction of the wave vector of the photon beam, fc, has to 
be specified. In the experiments, the single crystals were 
placed on the rotating holder so that the wave vector 
was parallel to the [101] crystallographic direction. Con- 
sequently, the spectra ((T|| -|- 2a±)/3 shown in Fig. (Sjare 
calculated with A;// [101] 4^ The two features of the exper- 
imental pre-edge are reproduced with or without the Is 
core-hole (see inset of Fig. ^ . The presence of the core- 
hole essentially provokes a 2.5 eV shift towards lower en- 




Energy (eV) 

FIG. 5: Influence of the core-hole effects in XANES calcula- 
tion at the Cr K edge in ruby. Inset: zoom of the pre-edge 
region. For pre-edge simulations, the calculated spectra in- 
clude both electric dipole and electric quadrupole transitions 
(the wave vector k is along the [101] direction of the trigonal 
cell of rubj4a). 

ergy of the pre-edge features, while the main rising edge 
is not shifted. A similar energy shift due to the core- 
hole was already observed for electric quadrupole tran- 
sitions at the Ti K pre-edge in rutile44 Unfortunately, 
here the shift induced by the presence of the core-hole 
is not large enough to give a perfect agreement with ex- 
periment. This problem may be due to the limit of the 
density functional theory within LDA to model the core 
hole-electron interaction for empty d states of 3d transi- 
tion metals. LDA+U calculations could possibly improve 
the agreement in the pre-edge region. Indeed in Ref.'*^, 
the authors show that LDA-I-U calculations performed on 
ruby essentially affect the positions of Cr d level relative 
to the valence and the conduction band. 



C. Comparison between experimental and 
calculated XANES spectra 

Fig. El compares experimental XANES data with the 
theoretical spectra obtained with the relaxed models. 
The (cr|| + 2iT^)/3 and — a± signals are displayed at 
the Ti K edge for the pink sapphire sample (FigEJi), at 
the Cr K edge for the ruby sample (FigEl)) and at the 
Fe K edge for the blue sapphire sample (FigO;). 

One observes that the experimental spectra present 
comparable overall shapes. The Fe K edge seems to be 
less resolved than the Ti and Cr K edges, probably due 
to the following two reasons: (i) the core-hole lifetime is 
shorter at the Fe K-edge (1.25 eV) than at the Cr K- 
edge (1.08 eV) and at the Ti K-edge (0.94 eV)j^ (ii) 
the instrumental resolution is estimated to 0.9 eV at the 
Fe K-edge while it is 0.6-0.7 eV at both the Cr and Ti 
K-edges. We also noticed similarities between the M K 
edges (M=Ti, Cr, Fe) in q;-A1203:M3+ and in a-M203 
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Energy (eV) Energy (eV) Energy (eV) 

FIG. 6: Comparison between experimental XANES data (dashed line) and calculated spectra obtained with the relaxed models 
(solid line): a) the Ti K edge in pink sapphire (Q;-Al203:Ti^''"); b) the Cr K edge in ruby (a-Al203:Cr^+); c) the Fe K edge 
in blue sapphire (a-Al203:Fe-Ti). An energy shift was added to valence band maximum in the calculations in order to match 
with experiments: 4966 eV at the Ti Kedge, 5991 eV at the Cr K edge and 7118 eV at the Fe K edge. 
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(see Ref4& for the Cr K edge in a-Cr203 and the Fe 
K edge in a-Fe203, and Refn^ for the Ti K edge in a- 
Ti203). At least to a certain extent, these spectra could 
be considered as the signature of a transition metal ion 
with a +III oxydation state in six-fold coordinated envi- 
ronment. However it is here necessary to perform calcu- 
lations to go farther in the XANES analysis. 

It can be seen in Fig. Elthat a good agreement is ob- 
tained between experimental data and calculated spectra 
carried out on the relaxed structural models. The agree- 
ment is especially good for ruby at the Cr K edge, where 
energy positions and relative intensities of the various fea- 
tures are well reproduced. It should be emphasized that 
not only the averaged (cr|| + 2ctj_)/3 spectra are correctly 
reproduced, but also the dichroic signal, from which more 
precise structural information could be extracted. Such 
a good agreement with angular dependent spectra has 
been also obtained with the same calculation method: at 
the Al K edge in corundum.^- and at the Si K edge in 
a-quartzi^ 

We have seen in Sec. IIII Al that the differences be- 
tween the relaxed and non-relaxed structures are essen- 
tially concerned with the interatomic distances within 
the coordination shell of the impurity. Therefore, it is 
interesting to test the sensitivity of XANES to such sub- 
tle structural modifications. Fig. [3 compares (T|| and cr^ 
experimental spectra with the corresponding calculated 
ones performed from both non-relaxed and relaxed mod- 
els. First of all, it is worth noting that the differences 
observed between calculated spectra are small but not 
negligible. For all compounds, peaks labeled g in <t\\ and 
labeled / and J in crj^ are too strong in the non-relaxed 
calculated spectra. Consequently, these features (all lo- 
cated at about 30 eV above the edge) are clearly corre- 
lated to the first interatomic distances around the impu- 
rity. The same conclusion can be made for the shoulders 
labeled A and B in cr^ and h in a\\ . Especially at the Cr 
K edge in ctx signal, one can see that the double-feature 
[D and E) is better reproduced by the relaxed calculation 
than by the non-relaxed one. 

In a previous study, the relaxed model was validated 
up to the second neighbors around the impurity by com- 
parison with EXAFS measurements. Here, thanks to the 
very good agreement obtained between experimental and 
relaxed calculated XANES spectra, the relaxed structural 
model is again validated and this time to greater extent. 



IV. CONCLUSIONS 

In this paper, we have shown that XANES can pro- 
vide precious information about the small structural re- 
laxation occuring within the atomic site of impurity in 
an aluminum oxide matrix. This has been achieved by 
(i) measuring angular dependent XANES spectra in or- 
der to probe the distortion of the impurity site, (ii) cal- 
culating structural model by ab initio energy minimiza- 
tion, (iii) calculating the XANES spectra from the previ- 
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FIG. 7: Comparison between (left) and a± (right) experi- 
mental data (thick solid line) and calculated spectra obtained 
with both relaxed model (thin solid line) and non-relaxed 
model (dashed line): a) the Ti K edge in pink sapphire, a- 
Al203:Ti^+; b) the Cr K edge in ruby, Q-Al203:Cr^+; c) the 
Fe K edge in blue sapphire, a-Al2 03:Fe-Ti. 



ous theoretical structural models using the full-potential 
pseudopotential self-consistent method of Ref The good 
agreement obtained between experimental and calculated 
spectra permits the validation of the structural models. 
We have also pointed out the importance of the elec- 
tronic potential construction to carry this study to a sue- 
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cessful conclusion. These results opens new applications 
in XANES analysis, in particular for dilute samples for 
which good quality EXAFS measurements are difficult to 
collect. 



Informatique Scientifique in Orsay, where all the calcu- 
lations of this study were carried out. 
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